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Abstract 

Arc object representation scheme which allows fast 
Euclidean distance calculation is presented. The object 
model extends the polytope model by representing ob- 
jects as the convex hull oj a finite set of spheres. An 
algorithm for calculating distances between objects is 
developed wh ich is linear in the total number of spheres 
specifying the two objects. 

1 Introduction 

The efficient calculation of distances between ob- 
jects is an important problem in the field of robotics. 
Many task planning and control algorithms require a 
knowledge of distances between objects. Collision de- 
tection and avoidance depends upon distance calcula- 
tion between a robot and other obstacles in the en- 
vironment. Many path planning algorithms use the 
shortest distance between objects, particularly those 
involving potential fields or iterative methods which 
test a path and then adjust it if there is a collision. 
All of these applications require an efficient distance 
computation method if they are to be used in a real 
time environment. 

In addition to applications in robotics, there are 
many computer graphics techniques which can benefit 
from distance or collision detection calculations. For 
example, collision detection can be used in ray- tracing 
to reduce the number of calculations required to ren- 
der a scene by only considering those objects which lie 
in the path of the rays. In a virtual reality simulation, 
where a person interacts with objects in the environ- 
ment, collision detection can be used to tell when the 
person is touching something. 

A distance calculation method that is often used 
involves checking for intersections between all of the 
faces, edges, and vertices of polygonal objects. This 
brute force method, although straightforward, is ex- 
tremely slow. Many different techniques have been 
proposed to speed up distance computation. Quad- 
tree representations can be used to quickly detect colli- 
sions between objects and a stationary set of obstacles, 
but become awkward and slow when the obstacles are 


moving. Some methods use extremely simple approxi- 
mations to all the objects such as calculating distances 
between bounding spheres. This is very fast, but quite 
inaccurate in most cases. Still others use norm one 
or norm infinity for the distance calculation, to avoid 
the computationally expensive square root operation. 
However, this approximation can give inaccurate or 
misleading results. 

A common method used to represent convex objects 
is the polytope model. In the polytope model, objects 
are represented as the convex hull of a set of points. 
Gilbert, Johnson, and Keerthi [1] present a powerful it- 
erative technique for calculating the distance between 
polytopes. The execution time of their algorithm is 
approximately linear with respect to the total num- 
ber of vertices in the two objects. The algorithm is 
roughly 0{n) where n is the total number of vertices 
in the two objects. Later. Gilbert and Foo [2] general- 
ize their result to include any convex object. However, 
in doing so, the advantages of having a single model 
for all objects is lost. 

This paper presents the concept ol a spherically ex- 
tended polytope, or s-tope. Instead of representing 
objects as the convex hull of a set of points, an s-tope 
is the convex hull of a set ot spheres. This allows a 
much wider range of objects to be easily described. Ob- 
jects with rounded surfaces can be represented with far 
fewer vertices than would be required with the poiy- 
tope model. Additionally, an efficient distance com- 
putation method based on that of Gilbert et al.. is 
developed and experimental results presented. The al- 
gorithm still exhibits the 0{n) execution time as does 
the poly tope method, but the spherical extension pro- 
vides a more general object model. The s-tope concept 
is a generalization of a spherical object geometry pre- 
sented by Tornero, Hamlin, and Kelley [3], [4] . 

2 Notational Conventions 

The notation used in this paper is similar to that 
used in [2], especially when dealing with concepts orig- 
inating in that or previous ([1] ,[oi.[6]) works. 

Vectors and points are represented by lowercase 
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boldfaced letters (x,y). and scalars by lowercase ital- 
ics, (a, b). Lengths of vectors and distances are defined 
as the Euclidean norm 

length of v = ||v|| = >/v - v, (1) 

where the ( •) is the inner, or dot, product. L nit vectors 
are written with a hat (v). 

Sets are denoted by uppercase boldfaced letters 
(P,S), and Co P is the convex hull of the set P. Other 
notation is introduced later in the presentation. 

To define the spherically extended polvtope model, 
it is first necessary to discuss the notion and notation 
of a polyope. A polytope is the convex hull of a finite 
set of points. When dealing with polytopes, the ver- 
tices which define the polytope are given the symbol 
p. Therefore, if P is a set of points {p 0 < Pi< ■ ; > Pnl* 
then the convex hull of P is an infinite set of points X 
expressed as 


X = 



: x = A t Pj p £ E P , Aq + + A rl 

i — 0 



3 S-topes 

An s-tope is the convex hull of a finite set of 
spheres. A sphere is denoted s = (c,r), where c 
is the center and r is the radius. A set of spheres 
S = {s 0 , s L , • • • , s n ) consists of the set of centers C = 
Coi c i, ■ • • , c n and the set of radii R = {^oTi . ■ ■ * > r n } - 
The convex hull of S contains an infinite set of spheres 
called the set of swept spheres which is expressed as 


CoS 


n 

(c, r) : c = ^AjCj, 

i=0 


where c* E C , r'i € R. 
and Aq + • • • 4- A rl = 1. 




*=o 


(4) 

(5) 


Note that the set of swept spheres does not include 
all possible spheres which can fit inside the s-tope, only 
those which are generated by Equation 3. An s-tope 
is represented in this paper by the symbol *5 where 
an s-tope formed from a set of spheres (s 0 , , • • ■ , 

is denoted as S 0 _ n . The spheres which define the s- 
tope are called the spherical- vertices of the s-tope. For 
simple s-topes (those made from four or fewer spheres) 
the spheres are listed explicitly. For example. Sijk is 
an s-tope with spherical-vertices s t , s : and s*. 

An s-tope is said to be overspecified if one or more of 
its spherical-vertices may be removed without chang- 
ing the convex hull. An s-tope which is not overspeci- 
fied is said to be a valid s-tope. 


4 The Simple S-topes 

As preparation for the discussion on distance calcu- 
lation, it is useful to have a detailed description of the 
s-topes formed from sets of one to four spheres. 

4.1 Sphere 

The simplest s-tope is a single sphere. An s-tope 
formed from a sphere, is denoted S, and is identical 
to that sphere. 



Figure 1: An example of a bi-sphere 


4.2 Bi-Sphere 

An s-tope formed from two spheres, and Sj , is 
denoted <S i; . A set of expressions describing the set of 
swept spheres which defines the bi-sphere is obtained 
by rewriting Equation 3 in the following manner. 


C,j(A t >) — Cj-rA^fCj—Cj), (bi 

r i';(A i; ) = r*» + “ r *)i 

where A^- E [0, 1]. (8) 


An example of a bi-sphere is shown in Figure 1. 
Hereafter, the conical surface generated between the 
two spheres is referred to as the side of the bi-sphere. 
There are several parameters which may be used to de- 
scribe a bi-sphere. The axis , v, ; , is the vector from c, 
to c j. The length of the bi-sphere is the magnitude of 
the axis vector. The convergence angle, a tJ of the bi- 
sphere is the angle between the axis and the side. The 
sine of a t j appears frequently and is given the symbol 
rjij. Another useful parameter is the direction along 
the axis. v :j . These parameters have the following re- 
lations, 


= c 7 - Ci, 


lij = 


H ~ r ‘ 

i v ui! 


tt, ; = sin ‘ f] tJ . 


(9) 

( 10 ) 

(H) 

(id) 


The vector n, is normal to the side ot a bi-sphere, 
projecting away from the surtace. and satisfies the fol- 
lowing constraint, 


n, ■ v.j = -rjij- (O) 

A bi-sphere is valid if and only if neither of the 
spheres s, or is completely contained within the 
other. Equivalently, a bi-sphere is not valid if \riij \ > 1- 
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Figure 2: An example of a tri-sphere 


4.3 Tri-Sphere 

A Tri-sphere is an s-tope formed from three spheres, 

, 3 j, and Sk> and is denoted S^k- Figure 2 shows an 
example of a tri-sphere. The set of equations describing 
all the swept-spheres in a tri-sphere is 

c ljk {\ ljk ) = Ci 4- Xijicj - Ci) + >,■*(€* - c,)(14) 
= r i + ( t'j — r,’) + — r i),(l*5) 

where A; ; E [0, 1], A,* E [0, 1]. (16) 

A tri-sphere has three bi-sphere edges, and two 
planar faces. The three edges are the bi-spheres S i; , 
$ik i and Sj k . 

The two planar faces are tangent to all three spheres 
3 ,*, sj , and Sky as well 35 to three edges *S, j , S t jt, 
and Therefore, the normals of the two faces must 
satisfy Equation 13 for each of the edges. That is. 


n 5 ■ v :; 

= -7y. 

(IT) 

n, v t k 

= -nuc, 

(IS) 

ri, * v ; - ; . 

= ~r)jk- 

(19) 


These equations imply an angle of 90 — ft»; between 
h s and vq Jt and an angle of 90 — a,-* between ri,, and 
Vjjt- Since the three equations are not independent, it 
is necessary to add the following constraint on ri, in 
order to calculate the two normals ri, and ri, , 

||ri, || = I- (20) 

Figure 3 shows the relationships expressed in the 
preceeding equations. 

A tri-sphere is not valid if one ot the spheres is com- 
pletely contained within the bi-sphere formed by the 
other two spheres. 

4.4 Quad-Sphere 

A quad-sphere is an s-tope formed from four 
spheres. The notation for a quad-sphere is the same 
as for the other objects. Namely, for a quad-sphere 
with spherical-vertices s t , s, , Sk> and s*. the symbol 
is Sijki ■ The equations for tfie quad-sphere are similar 
to those for the tri-sphere. 



Figure 3: Normals to the faces of a tri-sphere 


Cijki{Xijkt) = c j + A ij ( c ; “ + A ik( c k ~ Ci) 

+Aj/(c/ - Ci), 

r lJ k{X lJ k) = n + \ lj (r } - r t ) + A ifc (rjt - r t ) 
+Xil{n - r,*), 

where A t; E [0, 1], A,* E [0, 1], A u E [0,1]. (21) 

Just as a polytope with four or more points is gener- 
ally a polyhedral object with triangular facets, a quad- 
sphere (or an s-tope with more than four spherical- 
vertices) is an object with tri-sphere facets. Therefore, 
the normal to the surface of any of the facets may be 
calculated in the same manner as for a tri-sphere. 

A quad-sphere is not valid if any of its spherical- 
vertices are completely contained within the tri-sphere 
formed from the other three spheres. 

5 Polytope Distance Calculation 

The basic algorithm for calculating the distance be- 
tween s-topes is similar to an algorithm for polytopes 
developed by Gilbert, Johnson, and Keerthi [l], which 
in turn is based on work by Barr and Gilbert [o], [6] ; 
Since the s-tope distance algorithm is an extension oi 
the polvtope distance algorithm, an overview of the 
poly tope distance algorithm is presented here. For 
more details, refer to the works cited above. 

5.1 Polytope algorithm 

The poly tope distance algorithm takes advantage ot 
the fact that the shortest distance between two objects 
is the same as the shortest distance from the origin to 
the object formed by taking the Minkowski difference 
of the two objects. Therefore, given two sets of points 
Pi and Pt which define polytopes V L and 'Pa, a new 
polvtope P 3l can be formed from the set P 21 as 


P 21 = P 2 - P x = {p 2 - Pi : p L E Pi • P? € Pj} (22) 


The algorithm requires the ability to determine the 
point in a set which is farthest in a particular direction. 
The support function, h p, for a set of points P, is 

defined in [1] as 

/ip(v) = max{p v ; p G P}- (23) 

The point p which satisfies this equation is furthest in 
the direction v. 

The basic algorithm used is as follows. Given two 
polytopes defined by their sets of vertices Pi and Po, 
create a new' set of points P 12 using Equation 22. 

1. Start with an empty working set, P 

2. Guess at the direction of shortest distance, v. Use 
the direction between the centroids if you don t 
have any better information. 

3. Find the point p(v) in P l2 , furthest in the direc- 
tion v, and add it to the set P . 

4. Find the closest point to the origin, d, on the poly- 
tope formed from the set P using the polvtope 
distance sub-algorithm. 

5. Discard any points in P which are not necessary 
to define the point d 

6. Repeat steps 2-4 until |(^p l2 ( v ))" ~ ^ 

Tolerance 

Figure 4 show' four iterations of the algorithm on a 
planar polytope, 7^-6- 

For two polytopes V x and V 7t which have n and 
m vertices respectively, the polytope R l3 formed from 
the Minkowski difference of the two polytopes has n x 
m vertices. This means that for each iteration of the 
algorithm, n x m inner products must be evaluated 
to calculate the support function /ip 3 • However, the 
support function may also be calculated as 

h'p ( v) = /?-p a ( v) — h'p^ ( v). (24) 

This formulation only requires n-f ?? 1 evaluations of the 
inner product. 

Amazingly, the polytope distance algorithm almost 
alwavs converges within four or five iterations, even 
when there are hundreds of vertices. This makes the 
poly tope algorithm approximately 0(n + m) in com- 
putational complexity [1]. 

5.2 Polytope distance sub-algorithm 

The distance sub-algorithm [lj calculates the closest 
point to the origin on a polytope with one, two, three, 
or four vertices. Additionally, it returns the minimal 
number of vertices necessary to define this point. For 
example, if the closest point is one of the vertices, then 
that vertex is sufficient to define the closest point. It 
the closest point is on an edge, then two vertices (the 
end points of the edge) are necessary to define the clos- 
est point. Similarly, if the closest point is on a face of 
the polytope, three vertices are necessary to define the 
closest point. Note that the only time four vertices are 
required to define the closest point is when the point 
is actually inside the object. Therefore, after each it- 
eration of the distance algorithm, there are never more 
than three points in the working set. 



1. calculate shortest distance to initial working set 

(Pi) 

2. add p 4 to w'orking set and calculate shortest dis- 
tance to working set (p l5 p 4 ) 

3. discard p L and add p^ and calculate shortest dis- 
tance to working set (p 4 ,p 5 ) 

4. add p 3 and calculate shortest distance to working 
set (p 4 iPs iPs ) 

Figure 4: Example of the polytope distance algorithm 

6 S-tope Distance Calculation 

Many of the advantages of the polytope represen- 
tation are also present with the s-tope representation. 
For example, the Minkowski difference and the support 
function are both simply calculated. The Minkowski 
difference of the s-topes, S x and <S 3 defined by two sets 
of spheres Si and S 2 is the s-tope 5, 2 defined by the 
set of spheres Si 2 , where 


S 1 2 — {^12 — i 

( c 1 2 > r 1 2 ) : C 1 2 — c 2 — C i . 

(25) 

Ci E Ci 

. C’j G C 2 . 

(26) 

r l2 — r \ 


(27) 

ci € Ri 

. ro G R 2 } • 

(28) 


The support function for an s-tope is also very sim- 
ilar to that of a polytope. For an s-tope 5 defined by 
spherical vertices in the set S. the support function is 
computed as 

/ig(v) = max{(Cj • v) -f : (c, , G S}. (29) 



6.1 S-tope algorithm 

The distance calculation algorithm is very similar to 
the one presented above for poly topes. Given two s- 
topes defined by their sets of spherical-vertices S\ and 
S 2 3 create a new set of points Sfj using Equation 26. 

1. Start with an empty working set . S 

2. Guess at the direction of shortest distance, v. Use 
the direction between the centroids if you don't 
have any better information. 

3. Find the sphere s(v) in S l2 . furthest in the direc- 
tion v, and add it to the set S. 

4. Find the closest point to the origin, d, on the s- 
tope formed from the set S using the s-tope dis- 
tance sub-algorithm. 

5. Discard any spherical-vertices in S which are not 
necessary to define the point d. 



Figure 5 Distance from the origin to a bi-sphere 


6. Repeat steps 2-4 until |((fis [n ( v ))~ — (||d|[)'l < 
Tolerance 

6.2 S-tope distance sub-algorithm 

The s-tope distance sub-algorithm is actually a col- 
lection of four separate but interdependent algorithms. 
There are algorithms for finding the shortest distance 
from the origin to a sphere, from the origin to a bi- 
sphere, from the origin to a tri-sphere, and from the 
origin to a quad-sphere. The algorithms for calculating 
the distance to a sphere, bi-sphere, and tri-sphere, are 
based on work by Tornero, Hamlin, and Kelley [3], [4]. 

The sub-algorithm takes advantage of the fact that 
each of the collection is called under specific circum- 
stances by the s-tope distance algorithm. For exam- 
ple, the tri-sphere sub-algorithm will only be called 
after the bi-sphere sub-algorithm has been called. In 
general, if the distance is to be calculated from the 
origin to an s-tope. S n . with spherical- vertices s,-s ri . 
it can be assumed that the distance has already been 
calculated from the origin to the s-tope formed 

by the spherical-vertices and that all the 

spheres were necessary to define the closest 

sphere. Thus, it is useful to think of the sphere s n 
as being added to the s-tope and the distance 

computed to the resulting object. 

6.2.1 Distance to a Sphere 

In this paper, the shortest distance between the ori- 
gin and an object is written as 


6.2.2 Distance to a Bi-Sphere 

Calculation of the shortest distance between the ori- 
gin and a bi-sphere *5 £J is performed in several steps. 
First, the point on the axis of the bi-sphere closest to 
the origin is found by calculating , 



(31) 


Second, to satisfy the angular constraint, the value 
of A i; - is updated, 


= A, 


IMA iT>\\ 

ll v ‘; II 


tan a i; . 


Figure 5 illustrates these ideas. 


(32) 


Simple case (A^ £ [0, l]). For the simple case of 
A l; £ [0.1] the distance is calculated between the origin 
and Sij(Xij) using Equation 30. That is. 


VlST(Sij) = VIST (s t ;( A;j )) = I|c i; ( Ajj )|| — r ij (X ij ). 


Handling A out of range [0*1]. \\ hen the A is not 

in the range [0.1] . the shortest distance is calculated 
from the origin to one of the ends of S t , . In particular. 


VlSTiobject). 


X tJ >= 1 => VlST{5 tJ ) = V1ST(3j) 

X tJ <= 0 => VlS'T{S i: ) = VIST(3i) 


The distance between the origin and a sphere is cal- 
culated by the following formula, 

VlST{s) = ] | c J j - r. (30) 


However, since Sj is chosen by the distance algo- 
rithm and is closer (along the vector c t ) to the origin, 
A ij vvill never be less than zero, and that case does not 
need to be tested. 


o 




Figure 6: Pseudo-Polar Coordinate System 


6.2.3 Distance to a Tri-Sphere 

The process of calculating the shortest distance from 
the origin to a tri-sphere Sijk has three steps. First, 
the direction of shortest distance to the origin is calcu- 
lated. Then, the A’s (see Equation 14) are calculated 
which correspond to the closest sphere on the tri-sphere 
Sijk to the sphere The distance is then calculated 
from this sphere to the origin, unless the A’s are out 
of range (see Equation 16), in vvhich case further cal- 
culation is necessary. 


Calculation of surface normals. The direction of 
shortest distance, n,, when the A's are in range, is sim- 
ply the normal to one of the planar surfaces bounding 
the tri-sphere. These normals must satisfy Equations 
17 to 20. Two of the three independent constraint 
equations are repeated here. 


n ,‘Vjk = ~ r ljk , (35) 

n s Vji = -/?;/• (36) 

A pseudo-polar coordinate system is formed using 
v {; ,v t/l and h, : — v t; x v;.;. The vector n, is located 
in this system using two variables, 9 and o. 9 is the 
angle between n* and the plane formed by and 
Vfc/. 6 is the angle from v„ to the projection of n, 
in this plane, a represents tne angle between v, ; and 
v k{ . This coordinate system is shown in Figure 6. 

Using these three angles, the constraint equations 
may be rewritten as 

a - arccos(v,j • v,*). (37) 


COS 9 cos 0 — T)ij . 

(38) 

cos 9 cos (a + o) = i]i' K . 

(39) 

Solving for 0 and 9 gives 


(lik - r)ijtosa\ 

0 = arctan : , 

v TjijSina J 

(40) 

9 = arccos [ — ) . 

V cos o J 

(41) 


The two values of 9 specify the two faces of the tri- 
sphere. The vectors ii, are then constructed using the 
equations 


cos 8 sin(a — o) 
sin a 

(42) 

cos 9 sin 0 
sin a 

(43) 

g c = sin 9, 

(44) 

v^ + QikVik + 9cfi c . 

(45) 


Calculation of A’s. Let d be the length of the line 
segment joining the origin with the closest sphere in 
the tri-sphere S i; *. The following equation may be 
obtained from Equation 14, 


cl n 5 — Cjf A;j(Cj Cj) A|^(c^- c I( |. (46) 

Multiplying Equation 46 by n [; = i\ s xv tJ and n lk = 
h s x v { ^- gives the following two equations for A.; and 

a;;-, 


C»- 

Vi; n ik ' 


(47) 


Cj lljj 
v t r n i; 


(48) 


Shortest distance (simple case). Assuming both 
A’s are in the ranges specified by equation 16 these 
A's may be used in Equations 14 and 15 to calculate 
the closest sphere ^ ljk ( Kj Aik ) in the tri-sphere. The 
shortest distance between the origin and the tri-sphere 
5ij k is then calculated as, 

VlST{S ljk ) = VIST (a t ;*( A*; . A,-*;).) (49) 


6 



Table 1: Handling A’s out of range 


“ccr - 


A ij -r A t ;- 

VlS7(S iik ) = 

>= U 

>= 0 

<= i 

bimple case 

< 0 

>= 0 

<= i 

VXST(S ik ) 

>= 0 

< 0 

<= i 

VIST (Sij) 

>= 0 

>= 0 

> 1 

VlST(Sjk) 

< 0 

< 0 

<= 1 

\WS(VlST(Si } ), VXST{S lk )) 

< 0 

>= 0 

> 1 

MI N(2>XST(S U ), VlST(S jk )) 

>= 0 

< 0 

> 1 

ms (VIST (Si, ), VlST(S, k )) 


Handling A’s out of range. If the A’s obtained 
above fail to satisfy one or more of the inequality con- 
straints, then the shortest distance is between the ori- 
gin and one of the three edges of the tri-sphere S^k- 
Table 1 shows all the possible situations. In the table 
VXST(Sij) refers to the calculation of the distance be 
tween the origin and 
6 . 2 . 2 . 



Sij using the method of Section 


As was the case when calculating the distance to a 
bi-sphere, it is possible to eliminate some of the cases in 
Table 1. Sphere was added to Sij to lorm S^k, and 
that sphere is closer to the origin along the direction 
of shortest distance calculated for S tJ . Therefore, the 
closest sphere in *S tJ ^ will never be in Sij , and A,-* will 
never be less than zero. 

6.2.4 Distance to a Quad-Sphere 

The method for calculating the distance from the 
origin to a quad-sphere «S,y relies heavily on the fact 
that the distance has already been calculated to S, ; 
This insures that when the origin is projected along 
hj, the direction of shortest distance from the origin to 
Sijk , it will lie within the triangle formed by c, , c } , and 
Cfc. Adding the fourth spherical-vertex to the s-tope 
S : ji * creates three tri-sphere faces Sjji. Siu, and *S,w- 
The shortest distance to the s-tope o t jki is identical to 
the shortest distance to one of these three tri-spheres. 

To determine which tri-sphere face of the quad- 
sphere contains the closest sphere to the origin, the 
centers of all the spherical vertices of S^ki . and the 
origin, are projected into a plane perpendicular to n«, 
the result will appear something like Figure 7. The 
triangle in which the projected origin lies determines 
which tri-sphere face contains the closest sphere to the 


Figure 7: Projection of the centers of a Quad-Sphere 
into a plane 

Table 2: Determining the face which contains the clos- 
est sphere 


Xu 

Xil 

Xki 

Candidate tn-sphere 


< U 

V 

ji 

< 0 

7 

> 0 

$ikl 

? 

>0 

< 0 

£lh 


n> j = Vj/xh,. (54) 

n ki - v ki x n s (55) 

These vectors, define three planes which contain the 
point c; and are parallel to h, and to v t i . vji , and 
v*/ respectively. The dot product of the vector c ; with 
each of these vectors is taken, to give, 

Xi i = fiii'Ci, (56) 


origin. 

A simple procedure to determine which of the pro- 
jected triangles c ( -CjC{. c t c;c /. or c ; ctc/. contains the 
origin has been developed which avoids actually having 
to project the vertices into a plane. 

First, three vectors arc calculated as follows: 

if h c c t > 0 then 


I hi 

= n 5 

X v {/ . 

(50) 

n J< 

= i C 

X V;-/ . 

(51) 

nu 

= n, 

X V*,; 

(52) 

*ii/ 

= 

x n 5 . 

(53) 


Xji = fiji • ci . (57) 

Xki = (58) 

The tri-sphere containing the closest sphere is then 
determined from Table 2. 

The question marks in the table indicate that any 
value is allowed. The tri-sphere face chosen by this 
method contains the closest point to the origin unless 
the origin lies within the interior of the quad-sphere 
(or more specifically, within the interior of the tetra- 
hedron formed by the centers of the spherical-vertices 
of the quad-sphere). The origin is tested to determine 
whether it lies within the quad-sphere by calculating 
the normal to the plane of the centers of the candi- 
date tri-sphere, and then taking the dot product of a 
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with this normal. If the origin lies within the quad- 
sphere, the distance algorithm halts and a collision is 
reported. Otherwise, the distance to the quad-sphere 
is found by calculating the distance to the candidate 
tri-sphere using the method of Section 6.2.3. 

7 Experimental Results 

To determine the efficiency of the distance compu- 
tation algorithms, a program was written and tested 
using randomly generated objects with varying num- 
bers of vertices. 

Many of the equations used for calculating distances 
can be greatly simplified if all the radii of the spheres 
in a given object are assumed to be constant. In this 
case, the problem is reduced to calculating the distance 
between polytopes and then subtracting the sum of 
the radii of the two objects. The program was written 
to take advantage of this degenerate case, and runs 
faster when given sets of poly topes instead of s- topes. 
Consequently, two sets of experiments were performed. 
The first using all polytopes, and the second using all 
s- topes. 

In both cases, objects with from one to one hundred 
vertices were generated by randomly placing the ver- 
tices on the surface of a sphere. To generate polytopes, 
the vertices were all spheres of zero radius. To gener- 
ate s-topes, the vertices were spheres ot random radii. 
The times for calculating distances between pairs of 
these objects at varying separations was measured on 
a Sun Sparcstatiou 1. The results of these measure- 
ments are presented in Figures 8 and 9. In each figure 
calculation times are presented for the following cases: 
just touching, separated by a distance equal to the ap- 
proximate diameter of the objects, separated by three 
diameters, and separated by seven diameters. The two 
figures show the linear relationship between the total 
number of vertices and the calculation time. Addition- 
ally, it can be seen that the calculation time is greatly 
reduced when the objects are not in close proximity. 


8 Conclusions 

The s-tope model is a simple extension to the polv- 
tope model which allows many types of objects to be 
represented by fewer vertices. An efficient distance cal- 
culation method has been presented which has an ex- 
ecution time linearly related to the number of vertices 
in the two objects. The distance calculation time for 
s-topes is slightly greater than for poiytopes with an 
equal number of vertices. However, since many ob- 
jects may be represented with far fewer vertices using 
the s-tope model, there may be a net reduction in com- 
putation time. 
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